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Measurement of the Muon Anomalous Magnetic Moment to 0.7 ppm [jl[] 

Yannis K. Semertzidis®' 

^'Brookhaven National Laboratory, Upton, New York 11973 
for the g-2 collaboration [|| 

The experimental method together with the analysis method and results of the data taken in 2000 and prospects 
of the muon anomalous magnetic and electric dipole moment experiments are presented here. 


1. Introduction 

The g factor of a particle, defined as the ratio 
of the magnetic moment of the particle in units of 
its Bohr magneton, over the angular momentum 
of the same particle in units of h: 

magnetic moment 

eh/2mc /-. \ 

^ angular moment ’ ^ ' 

h 

is used to enhance our understanding of the un¬ 
derlying theory, like: 

1. It is used to indicate that the proton {gp = 
-1-5.586) and the neutron ( 5 „ = —3.826) are 
composite particles. 

2. The ratio gp/gn of -1.46 being close to the 
predicted -3/2 was the first success of the 
constituent quark model. 

3. The g-2 value of the electron is non-zero due 
to quantum field fluctuations. The agree¬ 
ment between the experimental and the the¬ 
oretical value is a triumph of both QED and 
of the experimental approach. 

4. The g-2 value of the muon is usually more 
sensitive to higher mass particles than the 
electron g-2 by the ratio of (m^/me)^ « 
40000. Therefore it is used widely to check 
the validity of the standard model and as a 
sensitive probe for physics beyond it. 

2. Theory 

The anomalous magnetic moment of the muon, 
defined as is the sum of QED, 


hadronic, and weak interaction contributions plus 
any new physics that may be present: 

a^(theo) = a^(QED) -|- a^(had) -I- a^(weak) + 
a^(new physics). (2) 

The theoretical values of the various contribu¬ 
tions, especially that of the hadronic one, are the 
subject of many papers |,|,|,|[0,|,§,@ : 

• a^(QED) = 11658 470.57(0.29) x 

10-1° (0.025ppm) (^. 

• a^(had) = 683.3(7.7) x 10-i° (0.66ppm), 

based on the e+e” value for the lowest- 
order hadronic correction of c^(had, 1) = 
684.7(7.0) X 10-1° (0.60 ppm) g since only 
these data can be directly related to 
a^(had) without further theoretical as¬ 
sumptions. The higher-order hadronic con¬ 
tribution of a;i(had, 2) = —10.0(0.6) x 

10-1° j^, and o^(had,lbl) = 8.6(3.2) x 
10-1° [Q. For completeness a;_j(had, 1) = 
701.9(6.2) X 10-1° (0.53 ppm) based on the 
r data g. 

• a^(weak) = 15.1(0.4)xl0-i° (0.03ppm) 

Then the standard model contribution 
a^(SM) = 11659169(7.7) x 10-i° with a rel¬ 
ative uncertainty of (0.66ppm). The Feynman 
diagrams of the second order weak contribu¬ 
tions are shown in Figure (|^). An example of 
a^(new physics) is the contribution due to SUSY, 
Figure (||), where the supersymmetric partners 
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of W and Z, the chargino and neutralino, are 
involved. Their contribution is estimated to be 

q^(SUSY) = 14 X 10"^° tan/3. (3) 

\ ^susy / 


3. Experimental Method 

A bunch of highly polarized muons with mo¬ 
mentum P « 3.09GeV/c is injected into a ring 
of 7.112 m radius, Figure (||), with 1.45 T dipole 
magnetic field of very high uniformity. The time 
distribution of the injected beam has an r.m.s. of 
25 ns. At « 90° from the injection point, the 
muon beam is kicked onto stable orbits by a fast 
magnetic pulse (kicker) Q. Within a super-cycle 
of 3.2 s there are 12 bunches, 33 ms apart from 
each other. Vertical focusing is provided by elec¬ 
trostatic quadrupoles |j^. Muons of that mo¬ 
mentum have a lifetime of « 64.4 /xs and are 
stored in the ring for ft! 1.4 ms. 


3.1. Principle of the g-2 Experiment 

The muon momentum vector in the lab frame, 
P, precesses under the influence of the electro¬ 
magnetic forces. The muon spin in its own rest 
frame, S, precesses under the influence of only the 
magnetic forces present in its rest frame. Muon 
decay violates parity in a maximal way; in the 
muon rest frame the most energetic electrons go 
along the muon spin direction. Then the energy 
of the electron is Lorentz boosted due to the muon 
momentum resulting in an energy which is mod¬ 
ulated according to the dot product of the two 
vectors 17 , S ■ P, which for a spin 1/2 par¬ 
ticle is an exact sine wave of angular frequency 
LOa, giving rise to the principle of g-2 frequency 
detection. 

The electrons, having on average less momen¬ 
tum than the stored muons, spiral inwards where 
they are collected by an inner ring of 24 electro¬ 
magnetic calorimeters [|^, and their signals are 
recorded by 400 MHz waveform digitizers (WFD). 
The acceptance of the detectors depends on the 
position of the muon beam and is estimated to be 
about 1/3 over all decays. In Figure we show 
the energy spectrum of the detected positrons 
when the muon spin in its rest frame is paral¬ 
lel to the muon momentum in the lab frame (the 



Figure 1. Second order weak contributions to 
involving the W, Z, and Higgs. The Higgs con¬ 
tribution is negligible for the present mass limits. 
The total relative weak contribution is estimated 
to be -1-1.3 ppm with an error of 0.03 ppm. 



Figure 2. SUSY contributions to involving 
the superpartners of W, and Z. Their contribu¬ 
tion, enhanced by the factor tan/3, may be larger 
than the weak. 



Figure 3. The g-2 muon storage ring is shown 
here with the thermal isolation cover (white) in 
place. 
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higher energy spectrum shown in red) and when 
the muon spin is anti-parallel (the lower energy 
spectrum shown in blue). Counting the number 
of detected positrons above an energy threshold 
of 2 GeV, one then gets a sine wave with the g-2 
frequency shown in the inset. 



Energy (GeV) 


Figure 4. The energy spectrum of the detected 
positrons in red (blue) for a muon spin paral¬ 
lel (anti-parallel) to its momentum in the lab 
frame. The red line corresponds to the higher 
energy positrons. The inset shows the number 
of positrons as a function of time for an energy 
threshold of 2 GeV. The red dots correspond to 
the highest number of detected positrons above 
2 GeV, whereas the blue to the lowest. 


The angular frequency Wa in the presence of 
both E and B fields is found 0 to be 


Wa =- 


ClnB -\- 


7 ^ — 1 


- a. 




( 4 ) 


where P-E = I3-B = Q is assumed. The electric 
field influences LOa because it advances the muon 
spin in its own rest frame (owing to Lorentz trans¬ 
formations, the electric field is partially trans¬ 
formed to a magnetic field) and the muon mo¬ 


mentum in the lab frame. In the case of a re¬ 
alistic detector with finite energy and time res¬ 
olution, there is some level of overlapping pulses 
(pileup), which produces a second harmonic of to a 
but also a first harmonic of tOa- More about the 
time development of these components and how 
we deal with them in our data is given later in 
the analysis section. 


3.2. Magic Momentum 

The muon anomalous magnetic moment has 
a value of approximately ^ which is 

the second-order (dominant) QED contribution. 
Therefore, for 7 « 29.3, the above equation (||) 
reduces to 


Wa = - B. 


( 5 ) 


In Figure (^ we show the spin vector getting 
ahead of the momentum vector as a function of 
time. 

The reason equation is valid comes from 
the fact that the g -2 precession is the difference 
between the muon spin precession in its own rest 
frame minus the momentum precession in the lab 
frame. The value of 7 « 29.3 corresponds to the 
case where the radial E field precesses the muon 
spin in its rest frame and momentum in the lab 
frame at the same rate. This is the reason for 
choosing a muon momentum of P « 3.09 GeV/c, 
a.k.a. “magic momentum” 1^, which for B = 
1.45 T corresponds to a radius of « 7.11 m. 


3.3. Origin of the Electric Field and Pitch 
Corrections 

Due to finite muon beam momentum width the 
cancellation is not exact and there is a need for 
a small electric field correction. Also, the muon 
momentum may not be exactly orthogonal to the 
external magnetic held, so the Lorentz transfor¬ 
mation of the lab electromagnetic helds into the 
rest frame helds are slightly modihed. This effect 
introduces a small “pitch” correction. Both the 
E held and pitch corrections are small, their sum 
being about - 1 - 0.8 ppm. 

3.4. Equation to Estimate 

A muon at rest and while in the presence of 
a magnetic held has its spin precessing with an 
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angular frequency given by 


2 \mu 


B. 


( 6 ) 


When combined with equation (||), this yields 


OJa 


R 


an = 


LOs - UJn 


X-R 


(7) 


Here i? = ^, and LUp is the angular frequency 
of the free proton in the B field, measured with 
NMR techniques |20 One then also needs to 


know the ratio A = —, the value of which is taken 


from another experiment. It follows that A = — 
with and /Xp the magnetic moments of muon 
and free proton respectively. 

The precision with which can be evaluated 
is determined by the accuracy of R and A. The 
value of A is determined by measuring the mi¬ 
crowave spectrum of the ground state of muo- 
nium finding A = 3.183 345 39(10). The pre¬ 
cision of R depends on the precision of tUa and 
ujp. The quantity coa is estimated by detecting 
the electrons produced by the decay of the muons. 
The statistical uncertainty of tOa is 


UJa LOa'yTpA^/Nl 

Here A is the g-2 oscillation asymmetry and N^, 
is the total number of detected positrons. 


special Fourier analysis technique is required to 
avoid introducing artificial effects into the width 
of the Fourier analyzed data |2^. The muon 
radial distribution so obtained is shown in Fig¬ 
ure i)- The muon momentum distribution is in¬ 
ferred from the radial distribution using the for¬ 
mula Pp = Pq X [1 + {R — i?o)(l — n)/Ro], where 
Pq is the central momentum, i?o = 7.112 m is the 
center of the muon storage region, R is the actual 
radial location of the muons, and n is the field 
focusing index. 

The fast-rotation structure is eliminated by 
randomizing the time of the injected muon beam 
with the average fast rotation period of ap¬ 
proximately 149.2 ns, the result shown in Fig¬ 
ure (|^). After the randomization the Fourier anal¬ 
ysis spectrum shows no structure at the fast ro¬ 
tation frequency. 

4.2. Coherent Betatron Oscillation Fre¬ 
quencies 

The muon storage ring lattice is shown in Fig¬ 
ure (§■ The muon beam is injected through the 
inflector [Q whose acceptance is smaller than 
that of the ring itself. This fact has as a result 
that the phase space of the betatron oscillations 
is not filled, resulting in betatron oscillations of 
the beam as a whole, called coherent betatron 
oscillations (CBO). Those oscillations are both 
horizontal and vertical and their amplitudes are 
described by: 


4. Beam Dynamics 


X = Xe+ Xo COS {uJxt + 


(9) 


4.1. Momentum Acceptance of the Muon 
g-2 Storage Ring 

The injected muon beam debunches due to mo¬ 
mentum dispersion with a lifetime of « 25 /iS. 
The time spectrum of the positrons detected 
by a single detector at early times after injec¬ 
tion is shown in Figure (^). The g-2 oscillation 
along with the slowly decaying fast rotation struc¬ 
ture is clearly shown. The rotation frequency of 
the stored muons depends on their momentum. 
Therefore their momentum distribution can be 
found by Fourier analysis of the arrival times of 
the detected positrons. The momentum accep¬ 
tance of the ring is narrow (0.6% total) and a 


y = yo cos {ujyt + 9y), (10) 

where Xe is the horizontal equilibrium radial po¬ 
sition away from the Rq = 7.112 m center of the 
ring, and xq (j/q) is the horizontal (vertical) CBO 
amplitude. The Wx = 2'iTfxi cOy = 27r/y, with 
fx ify) the horizontal (vertical) CBO frequencies 
given by 

fx ci /c(l - Vl - n), (11) 

and 

fy — 


( 12 ) 








5 



Figure 5. The muon spin S in its own rest 
frame gets ahead of its momentum vector P in 
the lab frame as time progresses. The detected 
energy spectrum of the electrons/positrons in the 
lab frame probes the product S ■ P. 


2000 Muon Radial Distribution 



Figure 7. The muon radial distribution obtained 
by Fourier analyzing the fast rotation structure of 
Figure (|). The central vertical line corresponds 
to the center of the muon storage region with 
Ro = 7.112 m and the two lines on either side are 
the edges of the storage region at Rq ± 4.5 cm. 



Figure 6. The bunched structure of the muon 
beam, clearly evident in the inset showing the 
8.5 — 13 /iS time range, distorts the g-2 oscillation 
at early times. The fast structure is decaying with 
a lifetime of « 25 /rs with only the g-2 oscillation 
remaining at later times. 



Figure 8. The detected positron time spectrum 
of Figure after randomizing the muon injec¬ 
tion time with a period of 149.2 ns; the fast rota¬ 
tion structure seen before, is eliminated. 
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where fc = =6.7 MHz, the cyclotron fre¬ 

quency, and n is the field focusing index. Equa¬ 
tions are not exact due to the discrete na¬ 

ture of the quadrupole coverage of the ring. There 
are small corrections given in reference [p^ . 


4.3. Field Focusing Index for the 2000 Run 

In Figure we show the vertical r'y ~ y/n 
versus the horizontal Vx — -s/l — n tune, along 
with the most important beam dynamics reso¬ 
nances of our weak focusing muon storage ring. 
The acceptance of a weak focusing ring has a 
rather wide maximum at n = 0.5, being at 90% 
of its maximum at n = 0.137. 

The n value is proportional to the voltage ap¬ 
plied to the quadrupole plates. Figure & Due 
to the presence of the magnetic field, and for rea¬ 
sonable residual pressures in the vacuum chamber 
of about 10“^ Torr, it would be very difficult to 
work at n = 0.5. This is so because there is a large 
number of low-energy trapped electrons circulat¬ 
ing in the quad region . A reasonable number 
to work with was around n = 0.136 — 0.137, in the 
middle between two relatively strong beam dy¬ 
namics resonances at n = 0.126, and n = 0.148, 
shown in Figure (10). The effect of the low-energy 
trapped electrons has been studied and is shown 
to contribute less than 0.01 ppm to the magnetic 
field. Their effect on the quadrupole electric field 
is negligible |l^ . 

From equation ( 0 ) we have that the horizontal 
frequency corresponds to fx — 466 kHz which is 
very near twice the g-2 frequency of « 229.1 kHz. 


5. Analysis of uja 

In the year 2000, we had a very successful run 
in terms of accumulating a lot of statistics. In 
Figure ( p^ we show the total number of positrons 
detected with E > 2 GeV as a function of time. 

The equation describing the ideal positron time 
spectrum is given by 

N{t) = No{E) [1+AiE) cos {cuat + ME))], (13) 

where A{E) corresponds to the g-2 oscillation 
asymmetry and (j)a{E) the g-2 phase, both of 
which depend on the energy threshold E] for 
E = 2 GeV, A = 0.4. 


superconducting 



Figure 9. The muon g-2 ring lattice indicating the 
four quadrupole regions, the kicker region and the 
inflector location. The muon storage region is a 
torus of radius Rq = 7112 mm and a cross section 
of a circle with a diameter of 90 mm. The inner 
ring of 24 detectors is shown as black squares. 



Figure 10. The tune plane with the most im¬ 
portant beam dynamics resonances of our weak 
focusing muon storage ring. The 2000 run was 
mostly taken at n ~ 0.137, in the middle of two 
relatively strong resonances, at n = 0.126 and 
n = 0.148. The 2001 run was taken at n = 0.142 
and n = 0.122. 
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Figure 11. A photograph of the electrostatic 
quadrupoles taken from the end of a vacuum 
chamber; the ring center is on the left. The ruler 
units are in inches. 


4 Billion Positrons with E > 2 GeV 



Figure 12. The number of positrons detected 
with E > 2 GeV as a function of time. The error 
bars are shown in blue and the points in red. 


The Fourier analysis of the residuals of the fits 
to the data using equation ([l^ ) is shown in Fig¬ 
ure dil)- The amplitude of the various peaks, 
especially that of /cbo = fx lias a large ampli¬ 
tude relative to the white noise present in the 
spectrum, implying that the CBO modulation is 
statistically very important. The two CBO side¬ 
bands are not of equal amplitude and in partic¬ 
ular not equal to ^AA^, with An the amplitude 
at fx- This precludes that Nq is the only CBO 
modulated parameter as was assumed in the 1999 
data analysis 

At early times the number of positrons shown 
in Figure ( p^ is of order 10 million per 149.2 ns 
bin. Therefore, very small beam dynamics ef¬ 
fects are important and noticeable in the least 
fits. Since the acceptance of the detectors 
depends on the position of the muon beam rel¬ 
ative to the detectors, the time and energy spec¬ 
tra of the detected positrons are modulated with 
the CBO frequency. As a result, the g-2 phase, 
asymmetry and the normalization Nq are all mod¬ 
ulated with the CBO frequency, thus becoming 
all time dependent. Since the CBO frequency 
is very close to twice the g-2 frequency, it turns 
out 1^ that the asymmetry and g-2 phase modu¬ 
lation are important effects that need special at¬ 
tention, consistent with M.C. simulations. The 
way they manifest themselves is by phase pulling 
the g-2 frequency with an oscillation period of 
To = (/^-/„)-/ ~ 130^s. The CBO modula¬ 
tion affects, as we said earlier, the energy spec¬ 
trum of the detected positrons and hence their 
average energy as a function of time. However, 
since the oscillation period is « 130 /iS, it was dif¬ 
ficult to distinguish it from other slow effects like 
gain change, muon losses, and pileup. In 2001, 
we took data at different n values, specifically at 
n = 0.142 and n = 0.122 corresponding to a hor¬ 
izontal CBO frequency of 491 kHz and 421 kHz, 
respectively. 

The time dependence of the CBO modulated 
effects is given by: 

1. No(t) = Nq [1-f An cos {2TTfxt + 0n)], 

with Tx the lifetime of the CBO modulation 
found from the data to be of the order of 
« 100 fjs. 










2 . A{t) = Ao[l + Aac cos (271 f^t + cj)A)], 

3. 4>a(t) = (j)ao + A^ cos (277 + (j)^) 


The time dependence of -/Vo(t) was found by 
strobing the energy spectrum of the 2001 
data at the g-2 frequency, see Figures ( pT|jl^ ). 
The fitting function used was of the form: 
Pi sin (27rt/P2 +-P 3 ) Next it was verified 

by looking at the residuals after the 5-parameter 
fit to the data, minimizing by fitting the data 
with various functions, M.C. simulations, etc. 
The time dependence used for A{t) and (j)(t) was 
not possible to verify from the data, only the M.C. 
simulations showed that they could not be too far 
off. 

The amplitudes of ^at, A a, and A^ are con¬ 
sistent with values from M.C. simulations. The 
values of the phases (j)N versus detector from the 
fits to the data are consistent with running from 
0 to 27r. That means that if the sum of all de¬ 
tectors is used, the amplitudes oi A^, A a, and 
A^ are reduced substantially, consistent with the 
values obtained with fits to the sum. 

The phase pulling of g-2 due to CBO is best 
depicted in Figure ( 0 ) where a straight line 
fit to fa = Wq/Stt versus detector gives fa = 
229 073.98 ± 0.14 Hz and a = 58.7/21, indi¬ 
cating that there is a consistency problem. A 
fit to a sine wave plus a constant, fa + Pi * 
sin (|^ X deC#-I-P 2 ), to the same data gives 
fa = 229 074.02 ± 0.14 Hz and a = 24.3/19 
with Pi = (1.20 ± 0.20) Hz, and P 2 = 1.91 ± 
0.17 rad. 

When the CBO modulation is included for 
the parameters No(t), A{t), and 4>a{t), the 
to a straight line fit of fa versus detector. Fig¬ 
ure ©, is = 23.9/21, and the average fa = 
229 073.92 ±0.14 Hz. 

5.1. Different Approaches to the CBO 
Modulation and Slowly Varying Ef¬ 
fects 

We have used several different approaches to 
analyze the data taken in 2000 : 

• Positrons with E > 2 GeV and a function 
including the modulation of No{t) and A{t) 
with fa:. 



Figure 13. Fourier analysis of the residuals to the 
fits to the data of Figure ( p^ using equation 0). 
fcBO = fx has a large amplitude implying that 
the CBO modulation is important. 



Time (}xs) 


Figure 14. By strobing the energy spectrum of 
the detected electrons of the 2001 data we were 
able to observe the time dependence of No{t). 
The high voltage was 21.7 kV, and n = 0.122. 
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Time (}xs) 


Figure 15. Same as in Figure ( p^ but for the 
case where the high voltage was 25.3 kV, and n = 
0.142. 



2 4 6 8 10 12 14 16 18 20 22 24 

Detector station 


Figure 16. A straight line fit to uJal^ir versus 
detector number when only the ideal 5 parameter 
function is used to fit the positron time spectrum. 
The for fh® straight line fit is = 58.7/21, 
and the average fa = 229 073.98 ± 0.14 Hz. A 
sine wave fit (see text) gives a good and a 
central fa = 229 074.02 ± 0.14 Hz, very close to 
the average. 


• Positrons in 200 MeV energy bins with 
1.4 < E < 3.2 GeV and a function includ¬ 
ing the modulation of A/)(t), A{t), and 
with fx- 

• Ratio method ; uja becomes independent 
of slow effects, e.g. muon losses. 

All the above methods gave results that are 
consistent among themselves within the expected 
statistical uncertainties due to the slightly differ¬ 
ent data used. 

In a side study, using positrons with E > 
2 GeV, we strobed the data at the horizontal GBO 
frequency fx making cua independent of the GBO 
parameters. Since fx is slightly higher than twice 
the g-2 frequency, one can recover all the infor¬ 
mation regarding the g-2 frequency, satisfying the 
Nyquist limit. Therefore the frequency, ampli¬ 
tude, and phase are recovered using mainly^] only 
the 5-parameter function and thus making no as¬ 
sumption as to the GBO functional form whatso¬ 
ever. This method gave, again, consistent results 
with the above methods. 

5.2. Pileup and other Systematic Errors 

When a positron arrives at the electromagnetic 
calorimeter and deposits an energy greater than 
approximately 1 GeV, it triggers the WFD con¬ 
nected to that particular detector. The WFDs 
are always running, and they were designed to 
keep in their memory more data, before and af¬ 
ter the pulse, than is necessary to reconstruct the 
positron signal that triggered them. This fact 
turned out to be of great help in dealing with 
the pileup pulses. Due to high rates in 2000, the 
overlapping of two positron pulses constitutes ap¬ 
proximately 0.5% of all detected pulses. We used 
the extra recorded data to reconstruct, on a sta¬ 
tistical basis, the time and energy spectrum of 
the pileup pulses which we then subtract from 
the data; see Figure 

The total systematic uncertainty in cUa is 
0.31 ppm §. The final/a = 229 074.11(14)(7) Hz 

^ There is also a slowly changing function multiplying the 
ideal function describing slowly changing effects, like muon 
losses, detector gain change with time, etc. Those are 
included to improve the overall but make no difference 
in the final fa value obtained from the fits. 
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Figure 17. A straight line fit to iOal^n ver¬ 
sus detector number when the ideal 5 parame¬ 
ter function including the CBO modulations of 
No{t), A{t), and is used to fit the positron 
time spectrum. The for the straight line fit is 
= 23.9/21, and the average fa = 229 073.92 ± 
0.14 Hz. 



_ ^^^ 

1 2 3 4 5 6 


E [GeV] 

Figure 18. The energy spectrum of the detected 
positrons with energy greater that 1 GeV at all 
times (thick line) and at only late times (thin line) 
when the rates are low. The dashed line shows 
the pileup-subtracted spectrum at all times. The 
inset shows the average energy of positrons be¬ 
fore (filled circles) and after (open circles) pileup 
subtraction for if > 2 GeV. 


(0.7 ppm) which includes a total correction of 
-|-0.76(^ ppm due to electric field and 
pitch [ [l6[|^ corrections. Those corrections have 
been studied in many different ways: analytically, 
particle tracking and spin tracking. For the lat¬ 
ter, the BMT 11^ equations were applied with our 
ring parameters, following the equations of chap¬ 
ter 11.11 (pages 556-560) of reference [0. The 
results from all the above methods are in agree¬ 
ment to ~ 0.01 ppm. 

6. Analysis of cOp 

The magnetic field of the muon storage re¬ 
gion was measured with an NMR trolley every 
two to three days while in between it was fol¬ 
lowed by 367 fixed NMR probes located on top 
and bottom of the vacuum chambers. The av¬ 
erage B-field, convoluted over the muon distri¬ 
bution, was obtained by two largely indepen¬ 
dent analyses of ujp ||]. The magnetic field mul¬ 
tipoles integrated over the azimuth of the ring 
for one trolley run out of 22 are shown in Fig¬ 
ure ([l^); the central field was 1.451274 T. The 
final ujp/ (27r) = 61 791 595(15) Hz. The total sys¬ 
tematic uncertainty in Up is 0.24 ppm |^]. 

7. Results 

In order to compute a^, both uja and tOp values 
are necessary. The analysis groups dealing with 
uJa and ujp worked separately and had applied se¬ 
cret offsets to their results until it was decided 
the analyses were finished. This avoided biases 
that could influence the choice of data selection, 
analysis method, etc. After the analyses seemed 
complete, during a collaboration meeting a secret 
ballot was held of whether or not we should reveal 
the offsets and compute . It was unanimous for 
revealing the offsets and compute a^, which we 
did. The results are 

R 

a + = - -- = 11 659204(7)(5) x 10"^°, (14) 

X — R 

with a relative error of 0.7 ppm. The experimen¬ 
tal world average becomes 


ap{exp) = 11659 203(8) x 10”^°, 


(15) 
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with again a relative error of 0.7 ppm. This ex¬ 
periment, like most of the muon experiments, is 
still statistics limited. 

In Figure ( ^ ) we give the recent BNL a^(ea;p) 
values, the average and the theoretical values 
based on the current standard model with the 
e“'"e“ data (solid horizontal line) and the t data 
(dashed line) for the hadronic contribution ||]. 

8. Discussion and Future Prospects 

The difference between the experimental value 
and the current theoretical prediction of is 

a;,(exp) - a^(SM) = 34(11) X 10"^° (16) 

which is a little over three standard deviations 
and it may indicate new physics. One should, 
however, wait for confirmation of the e+e“ data 
and understand the reason why the t based 
data imply a higher hadronic contribution than 
the e^e“ data, before making any claims as to 
whether or not new physics has been seen. The 
difference a^(exp) —a^(SM) = 17(11) x 10“^°, i.e. 
only 1.5 sigma, when the r data are used. On 
the experimental side, we have already accumu¬ 
lated about 3 billion electrons with E > 2 GeV, 
equivalent to a statistical power of approximately 
0.7 ppm, from our 2001 run with negative muons. 
We are currently analyzing those data and ex¬ 
pect to finish by early next year. We also have 
scientific approval for an extra four month pe¬ 
riod which the High Energy Division of DOE has 
not yet approved, though they should, in order to 
properly conclude the experiment. 

Assuming that the e“'"e“ data will hold, and 
if supersymmetry is responsible for the g-2 devi¬ 
ation, an EDM from similar quantum loops but 
with a phase giving rise to T and P violation is 
natural and likely. That would make the muon 
one of the best places to search for an EDM . 
Such an effort is currently underway |^. It 
promises to be sensitive to physics beyond the 
standard model |]^ , ^^2|j^ , |^ and to continue 
the exciting muon physics of the past and present. 
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